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'— 1 A search for non-resonant Z)^ — > 7r + /i + ^, and £)jj\ — > tt /i + fi + decays is performed 

I using proton-proton collision data, corresponding to an integrated luminosity of 1.0 

fb _1 , at -y/s = 7 TeV recorded by the LHCb experiment in 2011. No signals are 
observed and the 90% (95%) confidence level (CL) limits on the branching fractions 
are found to be 



in 
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B(D+^ 7T+/X+/0 < 4.1 (4.8) x 10~ 7 , 
2 B(D + ^ 7r~>"V + ) < 2.2(2.5) x 10~ 8 . 

> 7r-/i + /i+) < 1.2 (1.4) x 10~ 7 



X 



B(.D + ->- vr+/x + /tO < 7.3 (8.3) x 10" 



These limits are the most stringent to date. 
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1 Introduction 



Flavour-changing neutral current (FCNC) processes are rare within the Standard Model 
(SM) as they cannot occur at tree level. At the loop level, they are suppressed by the GIM 
mechanism fl] but are nevertheless well established in B + — > K + fi + fi~ and K + — > 
decays with branching fractions of the order 10 -7 and 10~ 8 , respectively (2|[3]. In contrast 
to the B meson system, where the very high mass of the top quark in the loop weakens the 
suppression, the GIM cancellation is almost exact in D meson decays leading to expected 
branching fractions for c — > u processes in the (1 — 3) x range [4|j6]. This 

suppression provides a unique opportunity to search for FCNC D meson decays and to 
probe the coupling of up-type quarks in electroweak processes, as illustrated in Fig. [I](a,b). 

The decay — > although not a FCNC process, proceeds via the weak 

annihilation diagram shown in Fie;. [1(c). This can be used to normalise a potential 
D + — > 7r + /i + yU~ signal where an analogous weak annihilation diagram proceeds, albeit 
suppressed by a factor |V^| 2 . Normalisation is needed in order to distinguish between 
FCNC and weak annihilation contributions. Note that, throughout this paper, the inclusion 
of conjugate processes is implied. 

Many extensions of the SM, such as Supersymmetric models with R-parity violation or 
models involving a fourth quark generation, introduce additional diagrams that a priori 
need not be suppressed in the same manner as the SM contributions pin. The most 
stringent limit published so far is B(D + 7r+/i + AO < 3.9 x 1(T 6 (90% CL) by the DO 
collaboration |8|. The FOCUS collaboration places the most stringent limit on the Df 
weak annihilation decay with B(D+ — > 7i + fi + fi~) < 2.6 x 10 -5 (9j. 

Lepton number violating (LNV) processes such as D + — > 7r~/x + /x + (shown in Fig. [1(d)) 
are forbidden in the SM, because they may only occur through lepton mixing facilitated by a 
non-SM particle such as a Majorana neutrino [l0] . The most stringent limits on the analysed 
decays at 90% CL are B(D+ -> n'^^) < 2 x 10" 6 and B(D+ s) < 1.4 x 10~ 5 

set by the BaBar collaboration [ll]. B meson decays set the most stringent limits on 
LNV decays in general, with B(B+ 7r-/z+/i + ) < 1.3 x 10" 8 at 95% CL set by the LHCb 
collaboration |l2] . 

This Letter presents the results of a search for Df, — > 7i + fi + {i~ and Df-. — > n~ 

decays using pp collision data, corresponding to an integrated luminosity of 1.0 fb _1 , at 
-y/s = 7 TeV recorded by the LHCb experiment. The signal channels are normalised to the 
control channels — > 7i + (p with <p — > fi + fi~, which have branching fraction products of 
B{D+ ->• tt+(0 = (1-60 ± 0.13) x KT 6 and 



B(D+ -> 7T+(0 fi+fi-)) = (1.29 ± 0.14) x 10~ 5 13 



2 The LHCb detector and trigger 

The LHCb detector |14| is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < r] < 5, designed for the study of particles containing b or c quarks. The 
detector includes a high precision tracking system consisting of a silicon-strip vertex 
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Figure 1: Feynman diagrams for (a,b) the FCNC decay D + — > 7r + ^ + /x , (c) the weak annihilation 
of a meson and (d) a possible LNV meson decay mediated by a Majorana neutrino. 



detector surrounding the pp interaction region, a large-area silicon-strip detector located 
upstream of a dipole magnet with a bending power of about 4 Tm, and three stations of 
silicon-strip detectors and straw drift tubes placed downstream. The combined tracking 
system has momentum (p) resolution Ap/p that varies from 0.4% at 5 GeV/c to 0.6% at 
100 GeV/c, and impact parameter (IP) resolution of 20 Lim for tracks with high transverse 
momentum (pt)- The IP is defined as the perpendicular distance between the path of a 
charged track and the primary pp interaction vertex (PV) of the event. Charged hadrons 
are identified using two ring-imaging Cherenkov detectors 15 . Photon, electron and 



hadron candidates are identified by a calorimeter system consisting of scintillating-pad and 
preshower detectors, an electromagnetic calorimeter and a hadronic calorimeter. Muons 
are identified by a system composed of alternating layers of iron and multiwire proportional 
chambers. The trigger 16 consists of a hardware stage, based on information from the 



calorimeter and muon systems, followed by a software stage that applies a full event 
reconstruction. It exploits the finite lifetime and relatively large mass of charm and beauty 
hadrons to distinguish heavy flavour decays from the dominant light quark processes. 

The hardware trigger selects muons with p^ exceeding 1.48 GeV/c, and dimuons whose 
product of pt values exceeds (1.3 GeV/c) 2 . In the software trigger, at least one of the final 
state muons is required to have p greater than 8 GeV/c, and an IP greater than 100 Lim. 
Alternatively, a dimuon trigger accepts candidates where both oppositely-charged muon 
candidates have good track quality, exceeding 0.5 GeV/c, and p exceeding 6 GeV/c. 
In a second stage of the software trigger, two algorithms select — > n + fi + fi~ and 
Dfs — > -n~ n + ji + candidates. A generic ji + fi~ trigger requires oppositely-charged muons 
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with summed px greater than 1.5 GeV/c and invariant mass, m(fi + fi~), greater than 
1 GeV/c 2 . A tailored trigger selects candidates with dimuon combinations of either charge 
and with no invariant mass requirement on the dimuon pair. 

Simulated signal events are used to evaluate efficiencies and to train the selection. For 
the signal simulation, pp collisions are generated using Pythia 6.4 17 with a specific 
LHCb configuration |18|. Decays of hadronic particles are described by EvtGen 19 . The 



interaction of the generated particles with the detector and its response are implemented 
using the Geant4 toolkit 20 as described in Ref. |21 



3 Candidate selection 

Candidate selection criteria are applied in order to maximise the significance of Df*. — > 
7T + /i + /i - and Dfs — >■ 7r~/! + /i + signals. The candidate is reconstructed from three 
charged tracks and is required to have a decay vertex of good quality and to have originated 
close to the PV by requiring that the IP x 2 is ies s than 30. The angle between the Dt^ 
candidate's momentum vector and the direction from the PV to the decay vertex, #d, is 
required to be less than 0.8°. The pion must have p exceeding 3000 MeV/c, px exceeding 
500 MeV/c, track fit % 2 /ndf less than 8 (where ndf is the number of degrees of freedom) 
and IP x 2 exceeding 4. Where IP y 2 is defined as the difference between the x 2 °f the PV 
reconstructed with and without the track under consideration. 



A boosted decision tree (BDT) [22] with the GradBoost algorithm 23 distinguishes 
between signal-like and background-like candidates. This multivariate analysis algorithm 
is trained using simulated D + — > 7i + fi + fi~ signal events and a background sample taken 
from sidebands around the D~K — > peaks in an independent data sample of 36 

pb _1 collected in 2010. These data are not used further in the analysis. The BDT uses 
the following variables: #d; X 2 °f both the decay vertex and flight distance of the Dt^ 
candidate; p and px of the candidate as well as of each of the three daughter tracks; 
IP x 2 °f the Djs candidate and the daughter particles; and the maximum distance of 
closest approach between all pairs of tracks in the candidate Df-. decay. 

Information from the rest of the event is also employed via an isolation variable, Ap T , 
that considers the imbalance of px of nearby tracks compared to that of the candidate 

MDU) - (Ep)t 

Pt " pt(z>5,) + (I>V (1) 

where pr(D^) is the px of the meson and (%2p)t is the transverse component of 
the vector sum momenta of all charged particles within a cone around the candidate, 
excluding the three signal tracks. The cone is defined by a circle of radius 1.5 in the plane 
of pseudorapidity and azimuthal angle, measured in radians around the Df g s candidate 
direction. The signal Df, decay tends to be more isolated with a greater px asymmetry 
than combinatorial background. 
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The trained BDT is then used to classify each candidate. An optimisation study 
is performed to choose the combined BDT and particle identification (PID) selection 
criteria that maximise the expected statistical significance assuming a branching fraction 
of 1 x 10~ 9 . The PID information is quantified as the difference in the log-likelihood 
under different particle mass hypotheses (DLL). The optimal cuts are found to be a BDT 
selection exceeding 0.9 and DLL M7r (the difference between the muon-pion hypotheses) 
exceeding 1 for each /x candidate. 

In addition, the pion candidate is required to have both DLL M7r and DLL^ less than 
and the two muon candidates must not share hits in the muon stations with each other 
or any other muon candidates. Remaining multiple candidates in an event are arbitrated 
by choosing the candidate with the smallest vertex \ 2 (needed in 0.1% of events). 

Candidates from the kinematically similar — > 7r + 7r + 7r~ decay form an important 
peaking background. A representative sample of this hadronic background is retained with 
a selection that is identical to that applied to the signal except for the requirement that two 
of the tracks have hits in the muon system. Since the yield of this background is sizeable, a 
1% prescale is applied. The remaining Df g , — > 7t + 7t + tt~ candidates are reconstructed under 
the Df* — > n + n + n~ and Df* — >■ 7r~/i + /i + hypotheses and define the probability density 
function (PDF) of this peaking background in the fit to the signal samples. 



4 Invariant mass fit 

The shapes and yields of the signal and background contributions are determined using a 
binned maximum likelihood fit to the invariant mass distributions of the Df\ — > 7r + /x + /i~ 
and Dt g s — > 7r~fi + [i + candidates in the range 1810 — 2040 MeV/c 2 . This range is chosen to 
fully contain the PDFs of the correctly identified D + and Df candidates as well as those 
of Df-, — > 7r + 7r + 7r~ decays misidentified as Dt, — > 7r + /x + /x~ or Df-, — > 7r~/i + /x + . 

The .D J-j — > 7i + fi + /j~ and — > 7r~/i + /x + signals are described by the function 

/(*) K 6X P ( o ixt'li ) ' W 

\2a 2 + (x -nYa LjR J 

which is a Gaussian-like peak of mean /x, width a and where «l(x < /x) and an(x > /x) 
parameterise the tails. The parameters of this shape are determined simultaneously across 
all bins (discussed below) of a given fit including the bin containing the control mode. 

The Dt, — > 7r + 7r + 7r~ peaking background data are also split into the predefined 
regions and fitted with Eq. |2j This provides a high-statistics, well-defined shape for this 
prominent background, which is imported into the corresponding subsample signal fit. 
The misidentification rate (the ratio of the yield in the signal data sample to that in the 
7r + 7r + 7r~ sample) is allowed to vary but is assumed to be constant across all bins in the fit. 
A systematic uncertainty is assigned to account for this assumption. 

A second-order polynomial function is used to describe the PDF of all other combina- 
torial or partially reconstructed backgrounds that vary smoothly across the fit range. The 
coefficients of the polynomial are permitted to vary independently in each bin. 
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The Df, — > 7i + ^i + fi~ and D~K — > n~ data are split into bins of m(/i + /i~) and 
m(7r~yU + ), respectively. The bins are chosen such that the resonances present in m(yU + /i~) 
in the case of Dfx — > -k + p + p~ are separate from the regions sensitive to the signal, which 
are in the ranges 250 - 525 MeV/c 2 and 1250 - 2000 MeV/c 2 . For the Dt a) 7r" / u + /i + 
search, the bins of m(n~fi + ) improve the statistical significance of any signal observed, 
as it is assumed that a Majorana neutrino would only appear in one subsample. The 
definitions of these subsamples are provided in Tables [T] and [2] Cross-feed between the 
bins is found to be negligible from simulation studies. 

The Dt g s —> 7r + n + \i~ and Dfs — > n~ ji + data are fitted independently, with the 
Dt g -. — > 7i + fi + fi~ sample being fitted in two parts due to the requirement of some of the 
software triggers that m(fi + fi~) exceeds 1.0 GeV/c 2 . A Df, — > fit excluding these 

trigger lines simultaneously fits the low-m( / u + / u _ ), rj, p/uj and bins. Another fit to the 
Df s -.—> -k + ijl + ijl^ data, including these trigger lines, is applied to the high-m (/i + /i~) and <fi 
bins. The <fi bin is needed as it provides a signal shape and normalises any signal yield. A 
simultaneous fit to the — > tt~ fi + p + data is done in all four m(ir~ p + ) bins. The <fi bin 
from the Dfs — > n + fi + fi~ data is again used to provide a signal shape and to normalise 
any signal yield. 

The invariant mass spectra together with the results are shown in Figs. [2] and |3j 
Background-subtracted m([x + fi~) distributions are obtained using the sPlot technique 124] 
and shown in Fig. |4 The signal yields are shown in Table jlj for — > decays, 
and in Table 2] for — >■ n~ decays. The statistical significances of the two 

observed peaks are found by performing the fit again with the background-only hypothesis. 
Significances of 6.1 and 6.2 a are found for D + — > tt + (i] — > and Df — > 7r + (77 — > 

decays, respectively. In comparison to Dfs^ n + ((f> — > /i + /i~), B(D + — > 7r + (r/ — > fi + fi~)) = 
(2.2 ± 0.6) x 10~ 8 and B{D+^ n + (r] = (6.8 ± 2.1) x 10~ 8 for the D + and D+ 

decays, respectively, and match those expected based on the D^ — > r]7r + and r\ — > p + fi~ 



branching fractions [13 . No significant excess of candidates is seen in any of the signal 
search windows. 



Table 1: Signal yields for the — > 7r + /i + /i fits. The 4> region yields differ due to the different 
trigger conditions. 

Trigger conditions 

Triggers without 
min+fx-) > 1.0 GeV/c 2 



W 

Bin description 
low-m(/i + ^~) 

V 

p/uj 



m(ii + n ) range [MeV/c 2 ] D + yield Df yield 

250- 525 - 3± 11 1± 6 

525 - 565 29 ± 7 22 ± 5 

565 - 850 96 ± 15 87 ±12 

850 - 1250 2745 ± 67 3855 ± 86 



All triggers 



high-m(/i + /i ) 



850 - 1250 
1250 - 2000 



3683 ± 90 
16 ± 16 



4857 ± 90 
- 17 ± 16 
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Table 2: Signal yields for the D+ } -»• 7r /^ + /U + fit. The ^ region from the — >■ ir + [i + n 
channel is used for normalisation. The particle 'x' is a n when referring to D^— >■ 7r~/i + ^ + data 
and a /U for — >■ ir + fi + [i~ data. 

Bin description m(fi + x~) range [MeV/c 2 ] D+ yield D+ yield 



850- 1250 2771 ± 65 3885 ± 85 

bin 1 250- 1140 7± 6 4± 4 

bin 2 1140- 1340 - 3± 6 3± 5 

bin 3 1340 - 1550 - 1 ± 6 6 ± 6 

bin 4 1540 - 2000 ± 4 4 ± 5 



5 Branching fraction determination 

The — > ir + fi + fi~ and D^— > branching fractions are calculated using 
B(£>+ w) = T7 ^ x -<4-J 1 x vr+(0 ^ V)), (3) 

where — > ir/i/i represents either — » ir + fi + fi~ or — >■ The relevant signal 

yield and efficiency are given by N D + _._.„„ and e D + ._.„„, respectively, and the relevant 

(s) ^ (s) W 

control mode yield and efficiency are given by N D + and e D + +{fj> ^ ll+li - ) , 

respectively. 

The efficiency of the signal decay mode and the control mode include the efficiencies 
of the geometrical acceptance of the detector, track reconstruction, muon identification, 
selection, and trigger. The accuracy with which the simulation reproduces the track 
reconstruction and identification is limited. For that reason, the corresponding efficiencies 
are also studied in real data. A tag and probe technique applied to B — > J/i/jX decays 
provides a large sample of unambiguous muons to determine the tracking and muon 
identification efficiencies. The pion identification is studied using D* + — > ir + (D° — > K~ir + ) 
decays. The efficiencies observed as a function of the particle momentum and pseudorapidity 
and of the track multiplicity in the event are used to correct the efficiencies determined 
by the simulation. The correction to the efficiency ratio is typically of the order of 2% 
in each m(/z + /i~) or m(-K~ n + ) region. Small relative corrections are expected since the 
signal and control modes share almost identical final states. 
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Figure 2: Invariant mass distributions for Df^ — > ir + p + p candidates in the five m(p + p ) 
bins. Shown are the (a) low-m(p + p~), (b) rj, (c) p/uj, (d) <p (including trigger lines with 
m(p + p~) > 1.0 GeV/c 2 ), and (e) high.-m(p + p~) regions. The data are shown as points (black) 
and the total PDF (dark blue line) is overlaid. The components of the fit are also shown: the 
signal (light green line), the peaking background (solid area) and the non-peaking background 
(dashed line). 
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Figure 3: Invariant mass distributions for — > ir in the four m(n /i + ) regions. Shown 

are (a) bin 1, (b) bin 2, (c) bin 3, and (d) bin 4. The data are shown as black points and the 
total PDF (dark blue line) is overlaid. The components of the fit are also shown: the signal (light 
green line), the peaking background (solid area) and the non-peaking background (dashed line). 
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Figure 4: Background-subtracted m(// + /x~) spectrum of (a) D + — > -K + n + fi~ and (b) Df — > 
ir + fi + fi~ candidates that pass the final selection. The inset shows the <j) contribution, and the 
main figure shows the r\ and the p/co contributions. The non-peaking structure of the low and 
high-m(/x + /x _ ) regions is also visible. 
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6 Systematic uncertainties 



Systematic uncertainties in the calculation of the signal branching fractions arise due to 
imperfect knowledge of the control mode branching fraction, the efficiency ratio, and the 
yield ratio. 

A systematic uncertainty of the order 10% accompanies the branching fraction of 
the control mode Df. — > 7r + (</> — > fi + and is the dominant source of the systematic 
uncertainty on the branching fraction measurement. 

A systematic uncertainty affecting the efficiency ratio is due to the geometrical ac- 
ceptance of the detector, which depends on the angular distributions of the final state 
particles, and thus on the decay model. By default, signal decays are simulated with a 
phase-space model. A conservative 1% uncertainty is determined by recalculating the 
acceptance assuming a flat m(fi + fi~) distribution. 

The uncertainties on the tracking and particle identification corrections also affect the 
efficiency ratio and involve statistical components due to the size of the data samples and 
systematic uncertainties inherent in the techniques employed to determine the corrections. 
The corrections depend upon the choice of control sample, the selection and trigger 
requirements applied to this sample, and the precise definition of the probe tracks. The 
binning used to weight the efficiency as a function of the momentum, pseudorapity and 
multiplicity is varied to evaluate the uncertainty. The uncertainty in the choice of phase 
space model is accounted for by comparing the efficiency corrections in the extreme bins 
of the m(/z + / u~) or m(n^fi + ) distributions. In total, the uncertainty due to particle 
reconstruction and identification is found to be 4.2% across all bins. 

Also affecting the efficiency ratio is the fact that the offline selection is not perfectly 
described by simulation. The systematic uncertainty is estimated by smearing track 
properties to reproduce the distributions observed in data, using Df, — >• 7r + (0 — > 
decays as a reference. The corresponding variation in the efficiency ratio indicates an 
uncertainty of 4%. Also, the trigger requirements imposed to select the signal are varied 
in order to test the imperfect simulation of the online reconstruction and 3% uncertainty 
is deduced. The sources of uncertainty discussed so far are given in Table [3j 

Final uncertainty on the efficiency ratio arises due to the finite size of the simulated sam- 
ples. It is calculated separately in each m(n + n~) and m(n~ n + ) bin. These contributions 
are included in the systematic uncertainties shown in Table |4j 

The systematic uncertainties affecting the yield ratio are taken into account when 
the branching fraction limits are calculated. The shapes of the signal peaks are assumed 
to be the same in all m(/z + /i~) and m(7i~~/i + ) bins. A 10% variation of the width of 
the Gaussian- like PDF, seen in simulation, is taken into account for variation across 
the bins. In each bin, the shape of the Df-. — > 7t + 7t + tt~ peaking background is taken 
from a simultaneous fit to a larger sample to which looser DLL^ criteria is applied. As 
simulation shows the shape of the PDF is altered by a DLL M7r requirement. A variation in 
the peaking background's fitted width equal to 20% is applied as a systematic uncertainty. 
The pion-to-muon misidentification rate is assumed to be the same in all bins. Simulation 
suggests that a systematic variation of 20% in this quantity is conservative. Contributions 
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to the yield ratio systematic uncertainty are found to increase the upper limit on the 
branching fraction by around 10%. 



Table 3: Relative systematic uncertainties averaged over all bins and decay modes for the control 
mode branching fraction and efficiency ratio. The number in parentheses refers to the Df decay. 



Source Uncertainty (%) 
Geometric acceptance 1.0 
Track reconstruction and PID 4.2 
Stripping and BDT efficiency 4.0 
Trigger efficiency 3.0 
B{Df s) -> 7T+ <f> (p + V~)) uncertainty 8.1 (10.9) 



Table 4: Total systematic uncertainty in each m(/i + ^ _ ) and m(-7r~/i + ) bin with the uncertainty on 
the control mode branching fraction, the efficiency ratio and the statistical uncertainty stemming 
from the size of the simulated samples added in quadrature. The numbers in parentheses refer to 
the Df decay. 

Bin description D^ s) -fr 7r + ^+^~ (%) D+ s) -fr tt~^ + (i + (%) 

low-m(/i + /i") n -8 (16.9) 

high-m(/i+/0 11.2 (15.5) 

bin 1 11.1 (17.0) 

bin 2 10.9 (16.4) 

bin 3 11.1 (16.0) 

bin 4 11.3 (16.0) 



7 Results 

The compatibility of the observed distribution of candidates with a signal plus background 



or background-only hypothesis is evaluated using the CL S method 25 , 26 . The method 
provides two estimators: CL S , a measure of the compatibility of the observed distribution 
with the signal hypothesis, and CLb, a measure of the compatibility with the background- 
only hypothesis. The systematic uncertainties are included in the CL S method using the 



techniques described in Ref. 25 , 26 



Upper limits on the D + — > 7T + // + /i~ and D + — > tx~ ji + jj, + branching fractions are 
determined using the observed distribution of CL S as a function of the branching fraction 
in each m(fi + ^) or m{ii~ bin. Total branching fractions are found using the same 
method and by considering the fraction of simulated signal candidates in each m(n + fi~) 
or m(ir~fi + ) bin. The simulated signal assumes a phase-space model for the non-resonant 
decays. The observed distribution of CL S as a function of the total branching fraction for 
D + — > is shown in Fig. [5j The upper limits at 90% and 95% CL and the p- values 

(1 — CLb) for the the background-only hypothesis are shown in Table |5j . 
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Figure 5: Observed (solid curve) and expected (dashed curve) CL S values as a function of 
B(D + — > ir + fi + n~). The green (yellow) shaded area contains the ±ler (±2<r) interval of possible 
results compatible with the expected value if only background is observed. The upper limits at 
the 90% (95%) CL are indicated by the dashed (solid) line. 



Table 5: Upper limits in each m(fi + fi ) and m(7r n + ) bin and total branching fractions at the 
90% and 95% CL and p-values for the background-only hypothesis. 





Decay 


Bin 




90%[xl0" 8 


95%[xl0- 8 ] 


p- value 






low-m(/i + /i _ 


) 


2.0 


2.5 


0.74 


D+- 




higli-m(/i + /i _ 


) 


2.6 


2.9 


0.42 






Total 




7.3 


8.3 


0.42 






low-m((U + /x _ 


) 


6.9 


7.7 


0.78 


Df- 


-> TT + fJ, + H~ 


high-m(/i + /i _ 


) 


16.0 


18.6 


0.41 






Total 




41.0 


47.7 


0.42 






bin 1 




1.4 


1.7 


0.32 






bin 2 




1.1 


1.3 


0.61 


D+- 


-> 7T~/! + /i + 


bin 3 




1.3 


1.5 


0.94 






bin 4 




1.3 


1.5 


0.97 






Total 




2.2 


2.5 


0.86 






bin 1 




6.2 


7.6 


0.34 






bin 2 




4.4 


5.3 


0.51 


Df- 




bin 3 




6.0 


7.3 


0.32 






bin 4 




7.5 


8.7 


0.41 






Total 




12.0 


14.1 


0.12 
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8 Conclusions 



A search for the — > 7i + fi + fi and — > ir fi + [i + decays has been conducted using 

proton-proton collision data, corresponding to an integrated luminosity of 1.0 fb -1 , at 
y^i = 7 TeV recorded by the LHCb experiment. Limits are set on branching fractions in 
several m(/x + /i~) and m(n^fi + ) bins and on the total branching fraction excluding the 
resonant contributions assuming a phase-space model. These results are the most stringent 
to date and represent an improvement by a factor of fifty compared to previous results. 
The observed data, away from resonant structures, is compatible with the background-only 
hypothesis, and no enhancement is observed. The 90% (95)% CL limits on the branching 
fractions are 

B(D + ^ 7r+/x + /0 < 7.3 (8.3) x 10~ 8 , 
B(Df-> ttVV") < 4.1 (4.8) x 10~ 7 , 
B(D + ^ tt>V + ) < 2.2 (2.5) x HT 8 , 
B(Df-> 7r-/x + /x + ) < 1.2(1.4) x 10" 7 . 
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